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Abstract 

If a leptoquark is produced at HERA as a narrow resonance, various ef- 
fects tend to broaden the measurable mass distribution considerably. These 
effects are discussed here, with special emphasis on initial- and final-state 
QCD radiation. A proper understanding is important to assess the signifi- 
cance of data and to devise strategies for better mass reconstruction. 
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Recently, the HI and ZEUS Collaborations at HERA have presented evidence for an 
excess of events at large Q 2 and x [|l], |2[. This could be nothing but a statistical fluke. 
Alternatively, it could be the first signal of the production of an s-channel resonance, a 
leptoquark (LQ). In order to test such a hypothesis, as further data accumulate, it is 
important to understand the production characteristics of an LQ at HERA. In particular, 
it should be noted that an expected narrow mass peak will be smeared by various physics 
and detector effects. In this letter we want to give a brief survey of the physics components 
that could contribute to this smearing, and estimate the magnitude of each of these effects. 
Previous studies of a similar kind have been performed with the Lego generator || , which 
is partly based on Pythia [|]]. We also want to give some examples of how to improve 
the reconstruction of the LQ mass. 

This is not a study on the physics implications of an LQ observation at HERA ||. 
Therefore we do not discuss the origin of an LQ in terms of an underlying theory, be that 
compositeness, supersymmetry or anything else, but stay with a purely phenomenological 
description of LQ properties |J. Even so, several options are possible, and we will restrict 
ourselves further. An LQ may have spin or 1, with differences in the decay angular 
distribution, but of little importance for the considerations in this letter, so we stay with 
the spin alternative. An LQ may have net fermion number (ql or qi) or ±2 (qi or ql) . 
Given that HERA has been running with an e + beam in recent years, the latter kind would 
be disfavoured by requiring a sea antiquark from the proton. It would have been favoured 
in the earlier HERA runs with e~, however, where the valence quark distributions could 
have been accessed. The ratio in parton distributions between valence and sea at x ~ 0.5 
being larger than the recorded luminosity ratio e + p/e~p, the previous non-observation |7| 
favours LQ's with vanishing fermion number. A leptogluon scenario is not excluded, but 
is not favoured ||. QCD consequences of a q£ LQ or a leptogluon are briefly mentioned 
later. 

Phenomenologically, the production cross section for a leptoquark in the process 
q + £ + — > LQ can be written as M 
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Here A gives the strength of the unknown Yukawa coupling; the alternative k parameter is 
normalized such that k = 1 corresponds to electromagnetic strength of the coupling. To 
first approximation x = MIq/s — corrections to this will be a main theme of the letter 
- and q(x, MIq) the parton distribution at a scale given by the LQ mass. 

The observed handful of candidate events per experiments, for an integrated luminosity 
of order 15 pb -1 and a detection efficiency around 50% (including a cut Q 2 > 15000 GeV 2 ), 
would suggest a production cross section of the order of 1 pb. For a ue + LQ of around 
200 GeV mass this corresponds to k ~ 0.01. The width of a scalar LQ, 

Tlq = ^M LQ = ^pM LQ , (2) 



then becomes only Tlq ~ 4 MeV. A de + LQ would be somewhat broader, but in either 
case the width is negligibly small on an experimental scale. Unless the total width is 
enhanced significantly by decays to exotic channels, the LQ is long-lived enough that it 
will form an LQ-hadron, made up of the LQ and an antiquark or a diquark. 

QCD radiation of gluons and QED radiation of photons in principle can occur in the 
initial state, off the LQ itself, and in the final state. Including all interference terms, a 
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complicated radiation pattern is then possible. The small width here offers a considerable 
simplification: interference terms and radiation off the LQ are suppressed for radiated 
energies above T L q || OT. The region E £ T L q ~ 10 MeV giving a very small contribu- 
tion to the total radiated energy, it is therefore sufficient to consider three main classes 
of corrections, (i) initial-state radiation, (ii) LQ-hadron formation, and (iii) final-state 
radiation. We will consider these effects one at a time, roughly in order of decreasing 
importance. 

The main effect of final-state radiation is that the hadronic system of the LQ decay 
acquires a mass in excess of the naive quark one and that, as a consequence, the energy 
of the recoiling lepton is reduced. The jet mass phenomenology is not very different from 
experience in hadronic e + e~ annihilation events. A standard parton-shower description 



TTJ tuned to LEP data predicts an average quark jet mass of about 30 GeV at 200 GeV 
energy. Taking over the same formalism for the LQ decays then gives an average mass 
(M q ) ~ 32 GeV, Fig. B], i.e. somewhat higher since the absence of a radiating final-state 
partner removes some phase-space competition. In a more detailed study one should 
include the explicit matrix-element information for the region of well-separated emission, 
which could introduce some modest dependence on the spin of the leptoquark. 
In the rest frame of the LQ the lepton takes an energy 

The rest-frame lepton scattering angle 9* is only very little affected by the QCD radiation 
|12|1 . (Furthermore, for a spin leptoquark, the inclusive decay distribution is isotropic in 
cos#* in any case.) Even though the standard DIS variables do not have the traditional 
meaning related to a spacelike boson propagator, they can be defined purely experimen- 
tally. The M q term above then propagates to give 

= ^(l-^), (4) 

_ <?_ _ M[ Q Q 2 
XBi ~ 2Pq ~ s Q 2 + M 2 • [b) 

Defining r to be the x value relevant for an LQ mass determination, 
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one finds that r > x^y In particular note that the denominator of the M q correction 
factor is Q 2 , not M 2 q. This means that low-Q 2 data would not be expected to show any 
peak in x By A typical cut Q 2 > 15000 GeV 2 reduces the (M q ) from 32 GeV to 29 GeV 
(eq. @))- The (M 2 ) « 1430 GeV 2 then gives an average correction factor in eq. @ of 
somewhat above 5%. If the experimentally reconstructed LQ mass is taken to be ^x-BjS, 
the original S function is smeared as shown in Fig. |2|. Since Q 2 = x^jl/s is an identity 
that follows from the definition of the respective variable, the alternative experimental 



measure yQ 2 /y gives the same smearing. 

Also photons may be radiated in the final state, both from the lepton and the quark. 
Owing to the smaller coupling, the amount of radiation is reduced compared with the QCD 
case above. Radiation almost collinear with the lepton occurs at a significant rate, but is 
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not resolved by the calorimetric definition of lepton energy and so is of little consequence. 
In our studies we choose not to resolve final-state photons below a 1 GeV invariant mass 
cut-off. 

We next turn to initial-state radiation, i.e. radiation off the incoming quark and 
lepton lines. This reduces the longitudinal momentum fraction carried by the reacting 
quark/lepton, and builds up a p±_ and a spacelike virtuality for it. 

Photon radiation off the lepton line is dominated by the almost collinear one. In this 
limit, one has that r = x q xe, where the X{ are the respective momentum fractions. In the 
leading-log approximation, the positron-inside-positron distribution is roughly Ul3 



D e Jx £ , M 2 ) « [3(1 - x.f- 1 with (3 = ^ f In ^ - 1 ] « 0.0575 . (7) 
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The eq. ([I]) is then modified to 
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o = J f dx q dxt q(x q , M[ Q ) Dt(x e , M L 2 Q ) 5(x q x e - r) . (8) 

One obtains (xg) ~ 0.99 or (x q ) ~ l.Olr. The steep fall-off of q(x q , M 2 q) dampens the 
tail to large 1 — x £ . Therefore typical experimental cuts on photon energy lost in the 
beam pipe do not make a big difference. The relation between the true and reconstructed 
masses is 

x 2 

r = x Bj - ~t~ — — -— , (9) 

1 - (l-x £ )/y 

so the mass shift may be in either direction, Fig. [| 

Radiation of photons and gluons in the initial state leads to a buildup of spacelike 
virtualities and transverse momenta for the incoming quark and lepton. The largest 
effects come on the quark side, Fig. [l|, where (Q q ) ~ 14.3 GeV and (p± q ) ~ 3.6 GeV, 
while (Qe) ~ 0.85 GeV and (p±e) ~ 0.15 GeV for the incoming positron, using the 
spacelike parton-shower formalism of ||14j| . (Numbers vary a bit depending on choice of 
parton distribution parametrizations etc.) The spacelike virtualities can be seen as a 
kinematical consequence of the p± kicks, so the Q q (i) and p± q (e) are strongly correlated. 
Note that (Q q ) M 2 : spacelike parton-shower evolution is constrained by the limited 
phase space and the steeply falling parton distributions at large x. That is, the probability 
that the daughter quark at xi comes from the branching of a mother quark at > Xi is 
related to the ratio q(xi_i) / q(xi) of parton distributions, and the integral of q(xi_i) over 
Xi-i > Xi is small. At each branching one has p 2 u ~ (1 — z^Q 2 , with Z{ = Xi/x^i, so 
from the argument above it is clear why also (p\ q ) (Q q )- 

The consequences of nonzero virtualities and transverse momenta can again be traced 
through the standard DIS variables. Also including the previous effects, i.e. M q ^ and 
xt 7^ 0, gives the relation 
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where a = 1 (10) 



Xi 

and b = x t Q 2 + M 2 + Q 2 q - 2p ±q k' L . 



(A correction formula in terms of another set of variables is found in The spacelike 
virtuality and p± kick may either increase or decrease the estimated LQ mass. In par- 
ticular, there is an explicit dependence on the azimuthal angle between the p± q of the 
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incoming quark (= p^lq i n this approximation) and the k' ± of the outgoing lepton. The 
resulting spread in reconstructed LQ mass is shown in Fig. [| 

A contribution is also given by the primordial kj_ of the parton-shower initiator 
inside the proton. This primordial kj_ is some combination of Fermi motion and initial- 
state radiation below the soft cut-off of the simulation program. A typical scale for the 
former would be (k\) ~ 0.25 GeV 2 , but even if this is increased to 1 GeV 2 the net effects 
are quite negligible, since the k± is added to the much larger p± of the perturbative parton 
shower. Only in some future study of the LQ p± distribution itself would the 1ow-p_llq 
tail be sensitive to the choice. 

The effects of LQ-hadron formation are not large, since the LQ is so massive that its 
motion is not significantly affected by hadronization effects. Data on b and c quarks are 
consistent with a non-perturbative fragmentation where the LQ-hadron retains an average 
fraction 

1 ~ 0.995 (11) 



M LQ + 1 GeV 



of the original LQ momentum, cf. |15] . The antiquark (diquark) contributes a typical 



fragmentation p± to the net p± of the LQ-hadron. The LQ has some Fermi motion inside 
the LQ-hadron, of the order of the antiquark constituent mass, that shifts the momentum 
of the decaying LQ. The net effects of all the hadronization contributions should be a 
more-or-less random momentum shift at or below the 1 GeV scale, i.e. negligible. A 
consequence of confinement is that the LQ mass cannot be defined unambiguously, so the 
detailed effects of hadronization have to be seen in the context of some specific scheme 
for relating the LQ-hadron mass to the LQ mass itself. 

The hadronic system will show some effects of the LQ lifetime. If an LQ-hadron is 
formed, the hadrons produced from the proton remnant and the initial-state-radiation 
partons decouple completely from those produced in the decay of the LQ-hadron. This 
means, e.g., that the charged multiplicity is independent of the LQ decay angle, i.e. of Q 2 
(apart from some small trigger-bias effects). By contrast, if the LQ is too short-lived for 
LQ-hadron formation, the initial- and final-state partons are connected into a common 
colour string . The string length and hence the multiplicity now depends on the angle 
between the proton remnant and the LQ hadronic decay products (just like in standard 
DIS processes, with W 2 oc Q 2 for fixed XBj). If the initial- and final-state shower effects 
are neglected, the variation at small Q 2 is significant, but is reduced to ~ 10% in the 
Q 2 > 15000 GeV 2 region. When showers are included, however, these tend to dominate 
the particle production characteristics in a global sense. A memory remains in the soft 
(low-momentum) region, e.g. defined by |p| < 1 GeV in the longitudinal rest frame of 
the LQ. (The lab frame is unsuitable, since the LQ motion in this frame introduces large 
spurious effects.) Over the standard high-Q 2 range, and with full inclusion of parton 
showers, the variation in charged multiplicity is here over 20%. In the unlikely event that 
the LQ is very short-lived, coherence effects will also appear in the parton-shower stage, 
and so the multiplicity variation will spread to all |p| ~ Tlq. It is therefore possible to 
conclude from the hadronic final state whether an LQ-hadron is being formed or not, 



given enough statistics. This situation is quite similar to the one for the top quark |L0 . 

All the physics aspects described above are included in the Pythia event generator 
@].[] The net distribution of reconstructed LQ masses is shown in Fig. |3], where all effects 
described above have been included. The bulk of the mass spread is given by the effects 



4 LQ-hadron formation is not part of the standard distribution. Note that MSTJ (50) =0 and MSTP (67) =0 
should be set to switch off interference between the initial- and final-state radiation. 



4 



in eq. flTOl). After correction with this formula, the remaining mass spread is the result 
of the k± kick and virtuality of the incoming lepton, the QED radiation off the outgoing 
lepton, and LQ-hadron formation. 

A q£ LQ would share most of the features described above. There is enhanced initial- 
state QCD radiation in those events where a branching chain q — > qg, g — > qq occurs. 
However, most events start with a q at the cut-off scale Qq ~ 1 GeV, and here gluon 
radiation off the q is suppressed by the steeply falling q(x) distribution. That is, parton- 
shower histories that allow a q to survive at large x and large Q 2 are biased towards 
lower activity (a bit more than those of the q), and so the (Qq) ~ 12.4 GeV rather than 
the 14.3 GeV of the q. For a leptogluon, an increased radiation in the final state follows 
directly from the higher colour charge, (M g ) ~ 45 GeV. In the initial state, there is again 
a suppression from the steeply falling gluon distribution. About half of the events contain 
a branching step q — ► qg, and this branching is biased towards higher virtualities since the 
quark line below the branching scale radiates less than the gluon above. The net result 
is a fairly large mean, (Q g ) m 31 GeV. A leptogluon would hence give a higher hadronic 
multiplicity than an LQ. 

So far we have concentrated exclusively on one method for LQ mass reconstruction, 
based on the measurement of the scattered positron (Q 2 , XBj and y). Alternative methods 
are used by the HI and ZEUS collaborations. However, they all share some common 
assumptions, such as massless partons and p± balance between the scattered positron 
and the recoiling hadronic system, and therefore do not recover the bulk of the smearing 
effects noted above. The optimal way to reconstruct a narrow(er) mass peak clearly is 
detector-dependent, and so must be worked out by the experimental collaborations. We 
here want to give some very general comments, however. 

With an ideal detector, and good separation between the LQ decay products and 
the hadrons from the proton beam remnant, the four-momentum of the LQ could be 
reconstructed to give the LQ mass directly. Such an approach suffers from the calori- 
metric smearing of hadronic momentum measurements. Since the parton-shower activity 
spreads out hadrons, there is also no clean separation of hadrons. Therefore, alternatively, 
one could start from the better-measured positron momentum and then add corrections 
according to eq. fllPf ), to give a mass estimate m reC on = \/ts. In practice, some of the 
required terms turn out to be rather difficult to reconstruct from the observable hadronic 
final state. The most feasible correction is for the final-state QCD radiation, which we 
have seen to be the largest individual souce of smearing, and which also tends to give a 
systematic bias in a mass determination. Two alternative methods have been compared 
for reconstructing M q . In the first, all particles with a (pseudo)rapidity below 1.7 (except 
the lepton) are used to calculate an invariant mass that is associated with M q . The cut 
is selected to minimize the misidentification of particles between the LQ-hadron and the 
remnant system. In the second method, a jet clustering algorithm is used to find jets 
with an E± > 10 GeV inside an R = yj (Ar/) 2 + (Aip) 2 < 1 cone, and an invariant mass is 
calculated for the particles belonging to the jet or jets. The jet threshold is here selected 
so that initial-state radiation is less likely to give rise to a jet. The reconstructed M q is 
then used in eq. (^) to give a corrected LQ mass estimate. The distributions in Fig. |]show 
that either correction factor does result in a narrower mass peak (the width is reduced 
from 11.1 to 8.2 GeV) and one that contains less systematic bias towards underestimating 
the correct mass (the average is shifted from 195.3 to 200.6 GeV). 

In summary, we have shown that several QCD effects are at play in the production and 
decay of a leptoquark. These effects are fairly well understood from our experience in other 
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areas of high-energy physics. It is therefore possible to describe in detail the production 
characteristics. This knowledge may be useful to obtain less biased and narrower mass 
peaks, and also to devise other tests that could help to distinguish between leptoquark 
production and ordinary DIS phenomena. 
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Figure 1: Distribution of M q (full), Q q (dashed) and pj_ q (crosses). 
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Figure 2: The reconstructed LQ mass distribution m rccon = y/XBjS for an input mass 
of 200 GeV, with a simple cut Q 2 > 15000 GeV 2 on the leptoquark sample. Results 
are shown with only one component active at a time: final-state QCD radiation 
(full), initial-state QED radiation off the positron (dashed) and initial-state QCD 
radiation (dot-dashed). 
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Figure 3: The reconstructed LQ mass distribution for an input mass of 200 GeV: all 
LQ events with m rccon = y/x^jS (dashed), Q 2 > 15000 GeV 2 with m recon = ^/xijs 
(full) and Q 2 > 15000 GeV 2 with m recon = yVi, eq. (0) (dot-dashed). 
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Figure 4: The reconstructed LQ mass distribution for an input mass of 200 GeV and 
with Q 2 > 15000 GeV 2 . Curves show mass based on xbj only (full), or corrected 
for M q either with a pseudorapidity separation (dashed) or a jet-finding strategy 
(dot-dashed). For details see text. 
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